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Chapter 10 Further investigation of integration 


We have two exceedingly important things to discuss in this chapter: integration over 
general subsets of R”, and change of variables. We first discuss a characterization of contented 
sets. 





A. Topological background 

We have often hinted at many of the following concepts. It is now time to make certain 
we understand them completely, and have useful names and notations for them. 

Recall the definition from Section 3A: if A is a subset of IR” and x9 € A, we say that 29 is 
an interior point of A if there exists r > 0 such that the ball B(xo,r) C A. 


S ~ 





DEFINITION. The set of all interior points of A is called the interior of A, and is denoted 
int A or int(A). 

As we have defined a set to be open if and only if all its points are interior points, we see 
that A is open => A = intA. 

Recall that we proved in Section 3A that the open ball B(x,r) is itself an open set. 


PROBLEM 10-1. Prove that every point of intA is in the interior of intA. That is, 
int(intA) = intA. 


PROBLEM 10-2. The preceding problem shows that intA is an open set. Prove that 
it is the largest open subset of A, in the sense that if B is an open set and B C A, then 
B Cc intA. 
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PROBLEM 10-3. Prove that 


int(AUB) D intA U intB; 


int((ANB) = intANintB. 


Prove that the inclusion expressed above may be a strict one. 


PROBLEM 10-4. Prove that the union and the intersection of two open sets are also 
open sets. 


The notion which is dual to interior will now be discussed. 





DEFINITION. Let A C R” and zg € R”. Then 2p is a closure point of A if for every 
0<r<oo, B(ao,r)N A is not empty. 


B(x,r) 


DEFINITION. The set of closure points of A is called the closure of A, and is denoted clA 
or cl(A). 


DEFINITION. A set is said to be closed if it contains all its closure points. 


PROBLEM 10-5. Prove that cl(clA) = clA. 


PROBLEM 10-6. Prove that the closed ball B(x,1) is a closed set. 


PROBLEM 10-7. Devise the analog of Problem 10-2 and prove the result. 
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PROBLEM 10-8. Devise the analog of Problem 10-3 and prove the results. 


PROBLEM 10-9. Prove that the union and the intersection of two closed sets are 
also closed sets. 


PROBLEM 10-10. Interior and closure are truly dual notions, for 


int(R” — A) = R”-clA. 





PROBLEM 10-11. Prove that A is open = > R” — A is closed. 
WARNING. Don’t misread this problem. A frequent error is to think that a set is open <=> 
it is not closed. 


We have discussed intA and clA. A third set will be of great interest to us: 


DEFINITION. For a given set A C R", the boundary of A is the set 


bdA = bd(A) 
clA — intA. 


PROBLEM 10-12. Let x € R” and A C R”. Prove that x» € bdA <=> for 


all 0 < r < o the ball B(zo,r) contains a point belonging to A and also a point not 
belonging to A. 





PROBLEM 10-13. Prove that bd(A) = bd(R” — A). 
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PROBLEM 10-14. Prove that R” is the disjoint union of the three sets 


intA, int(R"”-— A), bdA. 


PROBLEM 10-15. Prove that 


bdA = clANcl(R" — A). 


PROBLEM 10-16. Prove that 


A-—bdA = 
AUbdA = 


PROBLEM 10-17. Prove that 


bd(AUB) C bdAUbAB, 
bd(AN B) C bdAUbAB, 
bd(A—B) C bdAUbAB. 





We conclude this section with an important 
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LEMMA. Let A C R” and let IC R” be a special rectangle, and suppose that 
INA#@ and I-AFO. 


Then 
InbdAF 0. 





PROOF. The hypothesis implies there exist 
points x, y € J such that x € A and y ¢ A. 
Now the idea is to start at x and follow the 
line segment from x to y until we exit A y 
permanently: so we define the set 











S = {t|/0<t<1, (l—t)r+tye A}. 


Thus, S is a subset of the unit interval [0, 1]. Of course 0 € S and 1 ¢ S. We define tp = sup S. 
Then 0 < to < 1, and we set 


e Ift) =0, then z € A and the nearby points of [x, y| are not in A, proving z € A—intA C 
bdA. 


e If t = 1, then z ¢ A and there are nearby points of [x,y] which are in A, proving 
z€clA—AcC bdA. 


e If0 < to <1, then any neighborhood of z contains points of [x, y] which are in A (t < to) 
and which are not in A (t > to), proving z € clA — intA = bdA. 


In all cases, z € bd A. 
QED 





COROLLARY. If A is a proper subset of R” (i.e., A is not empty and A is not all of R"), 
then bdA is not empty. 





PROBLEM 10-18. Prove that in general a closed set is not equal to the closure of its 
interior. However, if A is the closure of an open set, then A = cl(intA). 


PROBLEM 10-19. Likewise, prove that if A is the interior of a closed set, then 
A = int(clA). 
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B. Topological characterization of contentedness 

Now we arrive at a tremendous result, showing that the property of a set’s being content can 
be expressed in terms of just the boundary of the set. The entire situation can be understood 
from the following equation. 


THEOREM. Let A C R” be a bounded set. Then 





vol(A) = vol(A) + vol(bdA). 


PROOF. This equation will be achieved by proving two inequalities. First, suppose P and 
Q are elementary polygons such that 






















































































PCAC, 
Then 

intP Cc intA CclA cclQ, 

and thus 

bdA Cc clQ — intP. 
Q 
P 

Therefore 

vol(bdA) < vol(clQ — intP) 


= vol(clQ) — vol(intP) 
= vol(Q) — vol(P), 
the last equality being due to the fact that the volume of a special rectangle is the same as 


the volume of its interior. It simply does not matter whether we use closed or open polygons. 
Now rearrange the inequality: 


vol(Q) > vol(P) + vol(bdA). 
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As P and Q are arbitrary, 
vol(A) > vol(A) + vol(bdA). 


The reverse inequality is apparently somewhat more difficult to establish. Suppose that 
R is an arbitrary special polygon which contains bdA. Choose a special rectangle such that 
RCT. Using the edges of R we construct a partition of J in such a way that R itself consists 

























































































of the union of certain special rectangles belonging to the partition. Thus we may express I 
as a nonoverlapping union of the form 


where each J; is aspecial rectangle. The crucial observation to make is that each open rectangle 
intJ; is either entirely contained in A or entirely contained in J — A. For if intJ; contained a 
point of A and also a point of J — A, the lemma of the preceding section would imply that 
int/; would contain a point of bdA, and thus of R. 
We conclude that the sets in the partition of J can be grouped into a nonoverlapping union 
as follows: 
l= - RUPUG, 


where P and Q are special polygons such that 
intP CA, inttQCI-—A. 
Thus 


vol(I) = vol(R) + vol(P) + vol(Q) 
= vol(R) + vol(intP) + vol(intQ) 
vol(.) + vol(A) + vol(I — A). 


IA 
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Since the theorem on p. 9-34 implies that 
vol(I) = vol(A) + vol(I — A), 


we conclude that 
vol(A) < vol(R) + vol(A). 


Finally, since R is an arbitrary special polygon containing bd A, 


vol(A) < vol(bdA) + vol(A). 


QED 
We have as an immediate corollary the definitive result about contentedness, 











THEOREM. A bounded subset of R” is contented <= its boundary has zero volume. 





PROBLEM 10-20. The theorem on p. 9-30 gives a very elementary proof that if A 
and B are contented, then so are AN B, AUB, and A — B. Give a different proof based 
on the theorem we have just proved. 


REMARK. We have now seen two proofs of the fact that if we deal with the collection 
of contented subsets of R”, then the operations of taking finite intersections, finite unions, 
and differences produce contented sets. This is described in set theory by saying that the 
contented sets form an algebra of subsets of R". The one finite set operation we can’t actually 


include is that of complementation, but this is only because we require contented sets to be 
bounded. 








C. Integration over more general sets 
Up to this point we have been integrating only over special rectangles. It is extremely easy 
to generalize to integration over arbitrary sets. 








DEFINITION. Suppose A is a bounded subset of R”, and A £, R is a bounded function 


defined on A. Let J be any special rectangle containing A, and define the new function I £.R 


by 
A= if EA, 





0 if cg A. 
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(The notation is intended to suggest that f has been extended to be zero outside of A.) Then 
f. is a bounded function on J, and we define 

1h 

Ce 


ae 
r= To 


(It is an easy matter to check that the definition of these new lower and upper integrals is 
independent of the choice of J; this same situation appeared also in Section 9H.) 
In particular, we can now simply write 


/ 1 = wld) 


[2 = vol(A). 


DEFINITION. We say f is integrable over A if these lower and upper integrals are the 
same, and we denote their common value as 


fe Lie T 


It is a straightforward task to use this definition effectively, so we will not have to say 
much about it. Here is an example, a generalization of the theorem on p. 9-33: 


THEOREM. Let A and B be bounded subsets of R", and let f be integrable over B. Then 





Lf 7 Id , fal 


PROOF. Let J be any special rectangle containing B and define f. to be the extension of f 
obtained by giving it the value zero outside B. Then f, is an integrable function on J. Then 
Problem 9-12 implies immediately that 


Ee = [ tetera f Ue felons) 


10 
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But we have the easy observations 
fl _ JF en BNA 
ena 0 outside BN A, 
f on B-A 
e et a 
Je~ Felana ' outside B— A. 
Thus we conclude respectively 
/ felBna = / ds 
— 4 “BNA 
[Ger fena) = ft 
i B-A 
QED 





PROBLEM 10-21. Suppose f is integrable over A and B is a contented subset of 
Prove that f is integrable over AN B. 








R”. 





D. Cavalieri’s principle 


This old and famous principle enables us “to find simply and rapidly the volumes of various 
geometric figures.” It’s an immediate consequence of Fubini’s theorem: 








PROBLEM 10-22 (Cavalieri). 
A(t) be the “slice” 


Let AC 














R” be contented, and for each t € R let 




















A(t) = {ye R"™ | (y, t) € A} 














(thus A(t) is a bounded subset of 


i 


(each of these integrals is of course 
empty for all sufficiently large |t]). 


(oe) 


vol,,(A) 


oe) 


vol,,_;(A(t))dé 


oe) 


VOln—1(A(t))dt 


R”-1). Prove that 


extended only over a finite interval of | 


oe) 





R, since A(t) is 
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PROBLEM 10-23. This is a favorite problem of basic calculus texts. Find the volume 
of the intersection of two right circular cylinders in R® which have the same radius r 
and whose axes are orthogonal. [You may assume the resulting set has the description 
ue +03 <r? 234232 <r] 








PROBLEM 10-24. Let B be a contented subset of R”~', regarded as the subset 
R" x {0} of R”. Let x € R®, with n™ coordinate z, =h > 0. Let C be the cone with 
base B and vertex x: 





























x 


C = {-t)yt+tz|0<t<1, ye Bh. 





It is a fact that C' is a contented subset of R”, but you need not prove that. Prove that 


vol, (C) = S96, at BA: 
n 











(The volume of a cone in R” is equal to 1/n times the altitude times the volume of the 
base.) 














Later in the present chapter we shall compute the volumes of balls B(x,7r) C R” for all n. 
We can actually already compute these, at least in principle. We introduce the following 








NOTATION. Ay, = vol,(B(0, 1)) 











It is an easy task to conclude that 


n 


vol, (B31) = aar™. 
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PROBLEM 10-25. Use Cavalieri’s principle to complete the following table: 


Q = 
ag = 
a3 = 
Aa = 
Os 


6 


PROBLEM 10-26. Calculate the volume in R? of the intersection of three cylinders: 
specifically, 


volg({z € R® |aj+az<1, ajt+a3<1, r3+23< 1}). 


PROBLEM 10-27. Calculate the 4-dimensional volume of the set 


{ER |aitaegtai<1, sitagtai<1, shtaptai <1}. 


PROBLEM 10-28. Calculate the 4-dimensional volume of the set 


{ore RY | ai toga <1, af bag ta, < 1}. 


PROBLEM 10-29. Calculate the 4-dimensional volume of the set 


fer |e 4a <1, aoa, <1). 
1 2 3 4 
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PROBLEM 10-30. Calculate the 4-dimensional volume of the set 


pe | a tat tai <1, 22+o7 <1). 
1 2 4 3 4 


PROBLEM 10-31. For any dimension n > 1 and any real number a > 0, let 
Sr(a) = {#@ ER" |0<a; for alli and 21 +---+2, < a}. 


Calculate vol, (S;,(a)). 


PROBLEM 10-32. A favorite calculus “stunt” is this problem: “A ball in R® has a 
cylindrical hole bored from it, the axis of the cylinder coinciding with a diameter of the 
ball. The set that remains has height h. What is its volume.” 

a. Use Cavalieri to calculate the volume. 

b. Why is this called a stunt? 


a 
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PROBLEM 10-33*. Generalize Problem 10-32 by assuming that the hole bored out 
comes from a right circular cone with axis coinciding with a diameter of the ball: 





Calculate the remaining volume as a function 
of the altitude h and the slant height 2. 


PROBLEM 10-34*. Generalize Problem 10-32 by assuming that the hole bored out 
is a paraboloid of revolution: 








Calculate the remaining volume as a function of h. 
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E. Elementary matrices 

We are now going to begin the final item of this chapter, the important change of variable 
formula for integrals on R". We shall first treat the case of linear variable transformations. 
This will be done in the following section, and in preparation for that we need to present an 
interesting feature of n x n matrices which we have not mentioned heretofore. 

As is our custom, when we think of a linear relation connecting xz and y € R”, we represent 
our points as column vectors and we write y = Tx, where T is of course an n x n real matrix. 
We certainly assume that T is invertible (det 7 4 0), as we are going to be interested in a 
one-to-one correspondence between x and y. 

We are particularly interested in two special kinds of what we shall call elementary matrices. 
In all cases below k and ¢ stand for fixed integers between 1 and n. 

Multiplying. This is any matrix of the following form. Let c 4 0. Let M be defined by 


























Mkk = C, 
Note that det M@=c and that M~ is also a multiplying matrix with c replaced by c™t. 
Mlustration (k = 1): 


c 0... O 
0 1 0 
00... 1 














Adding. In this case k 4 @ and let c € R. Let the matrix A be defined by 
ai = 1, 
ak = ¢, 

tie = 0 for all other (i, 7). 


Note that det_A = 1 and that A~ is another adding matrix, with c replaced by —c. Illustration 
(GL = 2): 


1c 0 0 
0 1 0 0 
000... 1 


Now notice the effect of multiplying a given matrix on the left by these two kinds of elementary 
matrices: 
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Multiplying. MT is the matrix obtained from T by multiplying row k by the number c. 
Adding. AT is the matrix obtained from T' by adding c times row I to row k. 


PROBLEM 10-35. Describe what happens to JT’ when T’ is multiplied on the right by 
an elementary matrix. 


It is quite amazing at first glance that just these two types of elementary matrices suffice 
to generate all invertible matrices, but the proof is easy: 


THEOREM. Every invertible matrix can be expressed as a product of elementary matrices. 


PROOF. Let T be an arbitrary n x n invertible matrix. We first want to maneuver to obtain 
the entry 1 in the (1,1) position. There are two possibilities. If t;; 4 0, then there exists a 
unique multiplying matrix M such that MT has the desired property. If t;,; = 0, then the fact 
that T is invertible implies that some t;; 4 0, and thus there exists an adding matrix A such 
that AT has the desired property. We describe both possibilities by saying that there exists 
an elementary matrix EF such that ET has the form 


1 
ET 


Then there exists a sequence of adding matrices A; such that 
1 
An—1An-2... A, ET = 
0 
By Problem 10-35, another sequence of adding matrices A, can be found such that 
1 0O...0 
fs 0 
An-1 eee A, ETA, aa. An-1 = 
0 
This new matrix we have manufactured contains an invertible (n — 1) x (n — 1) submatrix, 
to which we can apply the same procedure. This will not disturb the 1 or the 0’s which are 
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displayed above. Therefore, after n stages we arrive at a formula which can be expressed 
symbolically as 

E,Ey1... BoE \T F\Fo...PsiF, = I, 
where all the E’s and F’s are elementary matrices of the two kinds. Multiplying by their 
inverses in the correct way produces the formula 


P20 eb ee 





Since the inverses are also elementary matrices, the theorem is proved. 
QED 


Before presenting the integration application of this theorem about matrices, we pause to 
present an interesting proof of the profound multiplicative property of determinants. This 
fascinating proof is quite different from the proof we gave back in Section 3G. 


THEOREM. det(AB) = det A det B. 


PROOF. We are of course assuming A and B are n x n matrices. First we dispense with the 
easy singular case in which det A = 0. In this case the columns of A are linearly dependent 
(see p. 3-36). As the columns of AB are themselves linear combinations of the columns of 
A, the columns of AB are also linearly dependent. Therefore det(AB) = 0 and the desired 
equation holds. 

Now assume det A 4 0. Then we know that A can be expressed as a product of elementary 
matrices, 


A= Byes Eye 
But each of the two kinds of elementary matrices enjoys the elementary property 
det(EC) = det Edet C, 
as EC results from C' by a row operation as described above. Thus we conclude that 
det(AB) = det(E,E.... EB) 
= det EF, det(F.... E,B) 


= det E,...det E;, det B. 
Now use B = I to conclude also that 
det A = det Ey... Ex. 


Therefore, 
det(AB) = det Adet B. 


QED 
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PROBLEM 10-36. Prove that the representation produced in the theorem can be 
achieved with no more than n? elementary matrices. Why is n? “obviously” the least 
number of elementary matrices that can be employed in general? 


PROBLEM 10-37. Write out explicitly the representation of 2x2 matrices as products 
of elementary matrices. 


PROBLEM 10-38. Let T represent a 90° rotation: 


0 -1 
r= (02), 
Write T as a product of three adding matrices T = E, E23, and draw a sketch of the 


result of operating step by step on the square J = [0,1] x [0,1]. That is, sketch the 
parallelograms J, £3J, EoH3J, and TJ. 


F. Linear changes of variables 

In this section we are going to derive and prove a tremendous formula for changing dummy 
variables linearly in a Riemann integral over R”. The key ingredient in this formula will be 
the determinant of the matrix which encodes the linear change of variables. 














Before we launch into the actual statement and proof, let’s examine a nice special case. 
Let T be an invertible n x n real matrix, and express J’ in the usual manner of displaying its 
columns: 


T = (ty to te tn), 


where t; is of course the column vector 





Denote by J the special rectangle in R”, 
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Then we denote by T’J the set resulting from J by multiplying each member of J by T: 


TJ = {Tz|x2e J} 
— {PS | each O< a; < ; 
j=1 
= {Soa | each O< a; < i. 
j=l 
This set is therefore exactly the n-dimensional parallelogram with vertex 0 and edges t,,... , tn; 


see p. 8-1. We know from our intuztive treatment of volumes in Chapter 8 that 
vol(TJ) = |detT). 


That is, 
vol(TJ) = |detT| vol(J). 


We are now going to prove from integration theory that this formula holds in utter generality, 
in that J can be replaced by any contented set in R”. In fact, we shall prove a great deal 
more even than this. 

Before doing this, we clarify one small point. When one integrates on R! and makes a 
change of variable x = —y, one may write in the time-honored way 


[ seoa = - [1a 


This is not very convenient for us, as our usual condition a < b produces —a > —b in the 
“integral” on the right side. Thus we much prefer to write 


[ seoae = fey 


Notice the two things that occur with this viewpoint: (1) the change of variable y = —x 
transforms the interval [a,b] to the interval [—b, —a], and (2) the absolute value of dy/dx is 
the correct change of “differential”: 

















ic = ay, 


The above point of view is the convenient one for n-dimensional integration. In fact, it is 
convenient to regard the functions we are integrating as defined on all of R”, but equal to zero 
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outside some bounded set. Then we simply need not worry about the domain of integration, 
and we simply write the integral as 
/ feds. 
Rn 

Before starting the theorem, a lemma of a rather technical nature needs to be handled. It 
will help us in stating it if we employ a slight abuse of notation. If Tis an n x n matrix, let 
us also denote the associated function from R” to R” by the same letter T; that is, T’ is the 
function whose value at x € R” is Tx. (See Problem 2-77 for a discussion of the distinction 
that usually should be made.) As a result, we can write f oT for the composite funciton, the 
function whose value at x is f(T'x). 

You will see that virtually all the “hard” work involved in proving the desired result is 
contained in the following 














LEMMA. Let E be an elementary matriz, and let IR” 5,R" bea step function. Then fo E 


is integrable, and 

















Fly)dy = |det Bf f(B2)de. 
Rn Rn 
PROOF. Since step functions are linear combinations of indicator functions of special rect- 
angles, we may simply assume that f itself is such an indicator function: 
i. = ly where J = [a1,bi] x --- X [@n, Dn]. 
Note that (ly o E)(z) =1<— > Exe J == 26€ E"'J. Thus we have the simple formula 
ly of = lp-iy. 


We see then that the formula to be proved is 


| ly = jact 5) [ lp-1y; 
R” R” 


vol(J) = |det E|vol(E~*J). 


Of course we must also prove that 1,-1,7 is integrable, i.e., that E~'J is contented. 
In case E is a multiplying matrix, this is quite simple. As a typical case suppose 


in other words, 


ES ; : where a < 0. 
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Then : 
EVJ = 2.2] x [a2, be] x +++ x [an, dn], 
Ca 


so E~'J is itself also a special rectangle, thus contented, and 


vol(J) by — ay 
vol(E-1J) a ft 


a 


= -—a = |det E|. 





The case in which F is an adding matrix has a bit of a twist (really, a shear). As a typical 


case suppose 
1 


E = to 5 
Cc 1 


where the nonzero entry c is in the (n, 1) position. Then E~'J is not a rectangle, but is instead 
described by noting that x € E-'J => Exe J => 


IA 
IA 


ay Ly bi, 


Qn-1 Ln-1 ae 


= < 
Oy = Cry tay. Dy: 
The last condition can be rewritten in the form 
Gy, —CEy S Zn S ba — ch, 


and E~'J can be represented by a sketch in the x; — x, plane: 
b,-ca, 


ay,— Cay 
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First, this set is seen as the set contained between the graphs of two functions of x1,...,%n_1 
(namely, the affine functions a, — cx, and b, —cx,). As these functions are continuous, the set 
EJ is contented. This is an application of a slight modification of the theorem of Section 9I. 
Finally, Cavalieri’s principle gives the required formula for volume. The slice of E~'J for a 
given 7; =t, a; <t < }y, is the special rectangle 


[a2, by] Seok DK Caer Dn=i| x [an = ct, bn = etl. 
and the (n — 1)-dimensional volume of this slice is 
(be = a) nao (by F= An). 


Thus Cavalieri gives us 


vol(EJ) = / Wee te 


(61 — a1)(b2 — ag) +++ (On — Gn) 
= Vol J). 


Since det F = 1, the desired formula is proved. 
QED 


The hard work is over, and it is now a rather soft task to finish what we need to say about 
the linear case. 


THEOREM. Let T be a realn xn matrix with det T #0. Let R” 4, R be a bounded function 


which is zero outside some bounded set. Then 











f(y)dy = |detT|}  f(Tx)dz, 
R” R" 


and likewise for lower integrals. Moreover, f is integrable <> f oT is integrable, and in this 
case 


fly)dy = Jace 7) f (Laue 
R” Rn 


PROOF. The proof divides naturally into two stages. First we assume 7 = F is an 
elementary matrix. Then we consider any step function tT > f. Clearly, To E > fo EH. 
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As we know from the lemma, 7 0 F is integrable. Thus 

| doce) [f of < | det B| fo of 
| det E| / TOR 


/ t (also from the lemma). 


Since T > f is arbitrary, we obtain 


jot eifyon < fy ((¥)) 


That’s only an inequality. But we can employ the trick of replacing f by f o E~? in (x): 


| det aif Zi [i of. 


Another ruse: replace E by E~+ in this inequality, and multiply both sides by | det E], to 


obtain = = 
ie is jdct Bf fo B, 


This is the reverse of (*), so the theorem holds for the case of elementary matrices. 
The second stage of the proof is just the observation that if the theorem holds for two 
matrices 7 and 7), then it also holds for their product. For 


fi = jact | | 7 oT}. 


Then use the theorem for the matrix 75 and the function f o 7): 


[fen — jdet 7, [ f 07) oT, 


| det ff ie) (T,T9). 


Since the determinant is multiplicative, these two equations combine to produce 


iE = Jaot T35) ff 0 (1). 
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Finally, since any T’ is a product of elementary matrices, and the theorem is valid for elemen- 
tary matrices, the theorem is valid in general. 
QED 


PROBLEM 10-39. Prove that if A is any bounded set in R”, 


vol(TA) = |detT| vol(A), 


and likewise for lower volumes. Prove that A is contented <— > T'A is contented and in 


this case 
vol(T.A) = |det T| vol(A). 


PROBLEM 10-40. Prove that if ® € O(n), then 


vol(®A) = vol(A). 


PROBLEM 10-41. Let A be the ellipsoid 


where the semiaxes a@1,...,@, are any positive numbers. Prove that 


WOlL A) =: Oy Oy 1x5 Gy 


PROBLEM 10-42. Let M be a symmetric positive definite n x n matrix, and let A 
be the ellipsoid defined by 


Prove that 





EXAMPLE. Here’s a numerical example which nicely illustrates the power of these linear 
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changes of variables. Let A be the parallelogram in R? bounded by the four lines 


a and y = —~—4+2, 


Nog 
I 


Y= 22 and. y = 27—1. 


Compute (2, 3) 


/ / Pay, (4/3, 8/3) 
A 





We do not even need to draw a sketch 
of A, but here is one nevertheless: 
(—2/3, —4/3) 


Now define new dummy variables in a rather obvious way: 


uu = y—2z, 
v= Y=: 


This is a linear change of variables, represented as 


() = (i) Q): 


Notice that the determinant of the displayed matrix is —3/2. The four lines that form A are 
expressed in the u — v coordinates as 


= —1 and v = 2, 


u = 0 and u = —-1, respectively. 


Here’s the sketch in the u — v plane: 
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(-1, 2) (0, 2) 


(-1.-1) (0.-1) 


We know that dxdy should become | det T'|dudv, but T is really the inverse of the given matrix. 
Thus, dxdy becomes 2dudv. Finally, solving for x gives x = 2(v — u). Thus 


/ ih e’dxdy 
A 


3 


2 0 9 
jf ow —dudvu 
Lapis 3 
2 


0 
er’ dy . / e 24/3 day 
=i 


eer 


7 3 
(4/8 2p 2/8) (2/3 eit) 


4/3 _ 2/3 2/3 
(e € ) (e 1) ; 


WL wwl] rm wl bd 
No] w 


PROBLEM 10-43. Let A be the parallelogram in the x—y plane with vertices (1, —1), 


(4,0), (3,2), and (0,1). Compute 


[fy dady. 
A 


(Answer: 7/2) Can you now see how to find the area of A by inspection? 


PROBLEM 10-44. Let A be the triangular region in the x — y plane determined by 
the three lines « — 2y = 1, x+y =4, and y = 0. Evaluate the integral 





Ihy cau 
c= 2y 
A 


by using the coordinates u= 2+ y, v = x — 2y. 


(Answer: 17/9) 
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PROBLEM 10-45. Given n+ 1 points po, p1,...,Pn in R”, the simplex determined 
by them is the set 





j=0 j=0 


Thus a simplex in R! is an interval, a simplex in R? is a triangle, and a simplex in R? is 
a tetrahedron. Show that S can also be expressed in the form 


¢ = {rd — po) | all t; > 0 a <1}. 


j=l 
































Using Problem 10-24, show that 


1 
vol(S) = = | det(pi—Po P2— Po «* Pn Po) | - 


Each p; is written as a column vector. 
Pj 


PROBLEM 10-46. Show that the result of the preceding problem can also be ex- 
pressed as an (n+ 1) x (n + 1) determinant: 


voi(s) = [det (BB) |, 


nN: 





G. The general formula for change of variables 

The purpose of the present section is the understanding of the nonlinear generalization of 
the change of variables formula of the preceding section. This is one of the truly great results 
of calculus. It is the final tool which we need in order to make us adept at handling integrals 
in several variables. 

Important as it is, we choose to omit the proof. There are several reasons for this choice. 
First, knowing the proof will not help us at all in applying the result. Second, the context 
for providing an efficient proof is Lebesgue, not Riemann, integration. To be sure, we could 
give a proof using Riemann integration, but the technical details which would be required 
are actually far more involved than would be the case if we possessed the tool of Lebesgue 
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integration; the proof in the Lebesgue context is actually quite elegant. Third, the formula 
we are going to present is a very natural guess based on our understanding of the linear case, 
so much so that the result is quite believable. 

Thus our goal is to give a precise statement of the result, to provide several useful examples, 
and to provide exercises which illustrate the power of the result. We now give the exact 
hypothesis. 

Let A; and A», be bounded open contented subsets of IR”, and let 





Ay > As 


be a C! bijection whose inverse is also of class C?. 
We shall clearly see in some of the applications how these assumptions might be relaxed 
in various ways. 


THEOREM. Given the above situation, suppose 





A, 5R 











is integrable. Then f o® is also integrable, and 


f(y)dy = | f(®(x)) | det (2) | de. 
Ag Aj 
In this formula ®'(x) is the Jacobian matrix defined in Section 2H. 


REMARKS. The determinant of the Jacobian matrix of ® is often called the Jacobian 
determinant of ®. Notice that in the case of a linear function y = Tx, the Jacobian determi- 


nant is det 7. If we employ a sort of loose notation, regarding ® as representing y1,...,Yn as 
functions of 71,...,2%p, then the Jacobian determinant is precisely 
Oy. Oy 
Ox, — Orn 
det | : : 
Oyn. Oyn 
Ox, ae Oxn 


Notice also that in the linear case our result is precisely that of the preceding section. In 
that section the factor | det T| served as a type of magnification factor for volumes, as in the 
typical formula 

vol(['A) = |detT| vol(A) 


of Problem 10-39. 
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PROBLEM 10-47. Apply the theorem to the following situation: B C A, isa 
contented set and f is the indicator function on A, defined by 


fy) = 1 ©) €B. 


Show that the resulting formula is 


vol(®(B)) = [ leer) 


As a result of this problem we see that ® transforms subsets of A; to subsets of A» in 
such a manner that volumes of the transformed sets are calculated using integration over the 
initial set with a sort of local magnification factor | det ®’(a)|. You should be comfortable 
with this generalization of the linear case and should feel that approximating B with special 
polygons and using the continuity of | det ®’(x)| could reasonably be expected to provide an 
actual proof. 


P > 


aad @(B) 
A, 
A, 


POLAR COORDINATES. This is probably the most important of all examples. We use 
the usual formulas for polar coordinates in the x — y plane, 


2 = rcos?, 
y = rend, 
Here 0 < r < oo and @ is restricted to some interval of length 27. (Since r > 0 we are 


missing the origin, and if we think of 0 < 6 < 27 we are missing the positive x-axis. These are 
inconsequential details, as the missing sets have zero two-dimensional volume.) The Jacobian 
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determinant is 


cos@é —rsind 
det Ge 2) ~ 


so that when we change variables we replace drdy by rdrd0. (Notice the dimensional correct- 
ness.) Thus 


20 lee) 
/ j(@ yidedy = 7 | f(r cos 6,r sin 0)rdrdé. 
0 Jo 


EXAMPLE. We can now give a quick derivation of a recursion formula for a,,, the volume of 
the unit ball in R”. We do this by modifying Cavalieri’s principle to step down two dimensions 
at once. Thus 














vol, (B(0, R)) = a ; gee 
pee bey < 


/ / dxz3---dx,, | dx,dx 
x3 4-02 < R22? —a3 











x2 +a2<R? 
= // voln—2 (B00. \/R? — 2? — 3) dx,dxq 
x2+a2<R? 
n-2 
= | Qn—2(R? — 2} — 23) "2 dx dz, 
x2+a2<R?2 
olar coordinates n 
. I f On-2(R? —r an ~trdrdé 
r2)n/2 r=R 
= = 27An— os . 
—2 — aa 
20 
= —An- 2R”. 
nr 
Thus 
20 
An = —An-2 











We conclude therefore that 


2n, 27...,2m7 if n igs even 
= n n—2 4 ) 
Aan = 


an, on... 259 if n is odd. 
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In other words, 


ES 


An 


ols 


dl 


qn/2 

=a if n is even, 
n-1)/2 

=e! if n is odd. 


THE GAUSSIAN INTEGRAL. We calculated this integral once before, in Section 9G, as 
an illustration of Fubini’s theorem. We do it again here to illustrate polar coordinates. Again, 


A -{ ede. 
0 


as in Section 9G, define the number 


Then 


A? 


Fubini 


polar coordinates 





Thus A = 5. This was our second computation of the Gaussian integral. We'll see a third 


in Section H. 
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SPHERICAL COORDINATES FOR R?°. As we discussed in Section 6D, we use 





x = rsinycos8, 
y = rsingsing, 
2 = PWOsSe. 
Here 0 < r < =, 0< yp <7, and 0 < @ < 27. The Jacobian determinant is (if we use the 
ordering r, y, 0) 
sinycos@ rcosycosé —rsinysiné 


det | sinysin@ rcosysin@ rsinycosé 
cos yp —rsiny 0 


This can be calculated directly to be r?siny. An alternate point of view is that the columns 
of the matrix are orthogonal vectors whose lengths are, respectively, 1, r, and rsiny. The 
determinant is + the volume of the rectangle with these three sides, and is thus +r? siny. We 
don’t care about the sign in the context of integration, as we use the absolute value. Notice 
the dimensional correctness of the formula 


drdydz = r’sin pdrdpdd. 


As an example we compute again the volume of the unit ball, 


a3 = / dadydz 
x2 +y2+22<1 
QT T 1 
= | | | r? sin ypdrdpdé 
0 o Jo 
Qn T 1 
| ao. | sinyay: [ r’dr 
0 0 0 
1 


— O97 -Q-2 
3 


Ar /3. 











PROBLEM 10-48. The formula for spherical coordinates in R* is given in Problem 6— 
16. Show that 





dvoly = r°sin? (2 sin yi drdy dy2dé. 


Use this formula to recalculate a4. 
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INTEGRATION OF SPHERICALLY SYMMETRIC FUNCTIONS. Our experience 
with R?, R?, and R* reveals a definite pattern. We could write down spherical coordinates for 
IR” in terms of r and n—1 angles; these angles are y1,..., Yn—2, and @. Here each 0 < y; < 7, 
while 0 < 6 < 27. These must result in a formula of the form 


























dvol, = r”lo(p1, +> Yn-2)drdy, ...dYy_2d6, 


where we could calculate g if we wished. 











Now suppose IR” 4, Ris spherically symmetric with respect to the origin. This means 
that f(a) depends only on ||z||. By abuse of notation we choose to write f(||z||). Then the 
integration formula becomes 


love) T wT 27 
/ F(\le\)dri---dz, = | | af | f(r)r™"gdrdy,...dpp_2dO. 
R” 0 0 0 0 


In this integration we can use Fubini’s theorem to split the r integration and the angle inte- 
grations. The result is 


dx4+-*din = C * Pye an, 
[ flleldar--de, = ¢ f° ser) 


where C is a certain constant. 


PROBLEM 10-49. Choose a certain example for f in the above formula to conclude 
that C = na,,. The result is therefore 





[Heide ton = nam f ” p(n) Nar, 
R" 0 











PROBLEM 10-50. The unit ball B(0,1) can be inscribed in the cube [—1, 1] x --- x 
[—1,1]. The ratio of the volume of the ball to that of the cube is 2~-"a,. Prove that 
this ratio is a strictly decreasing function of n. Also prove that the ratio tends to zero as 
n — oo. In fact, prove that lim a, = 0. 


n— Oo 


H. The gamma function 
This section is not actually about changing variables in integrals, but instead is a natural 
conclusion to what we have learned about the number a,,, the volume of the unit ball in R”. It 
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is somewhat clumsy to have to use two different sorts of expressions for a, depending on the 
parity of n, as we have done in the preceding section. We shall soon derive a single formula 
that works for all n. 

The key lies in the introduction of a function which elegantly interpolates the factorial 
function, which is itself defined only for nonnegative integers. This function is rather artificially 
called the gamma function. The Greek letter gamma has no real significance in this context, 
just as the Greek letter pi in itself has no particular relation to circles. 





DEFINITION. The gamma function is the function (0, 00) +R given by the expression 
Cia) = i oe "di, 0a =< oo: 
0 
We remark that 0 < ['(a) < co. The defining integral is “improper” at infinity, but the 


exponential decay gives convergence; also the integral is improper at zero if 0 < a < 1, but 
comparison with ¢*~! gives convergence. 


PROBLEM 10-51. Use an integration by parts to prove that 








T(ia+1) = al(a). 








It is clear that P'(1) = 1, so that Problem 10-51 implies that 











I(a+1) = a! al a ee 





Thus we can assert that ['(a), defined for 0 < a < o, interpolates the factorial function 
(a — 1)!, defined only for a = 1, 2, 3,... . 


PROBLEM 10-52. Use a change of variable to show that 
ria) = 2 | we de, 
0 


Conclude that 














The equation for T (4) is truly fascinating! Just think: a rather natural interpolation of 
the factorial function produces such an interesting outcome for the “factorial” of —1/2. 
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Now look at the formula for a,,, as given on p. 10-31. If n is even, the denominator is 





Sr 
3 
N] | 
Nw 
rae 
| 
fe 
Nm] s 
n4/ 
| 
ie 
aS 
Nm] Ss 
+ 
rae 
nn" 


If n is odd, it is instead 





N | 
3 
Ww] | 
ie) 
ON es 
II 
ej 
GX 
w|s 
+4 
— 
NY 
TS. 
as | 
fs 
Nl rR 
4 


Thus in all cases we have the single formula 





gn/2 


r(24+1) 


An = 











Incidentally, some books write a! = [(a +1) for all —1 < a < oo. With such a definition, 
then we have the elegant expression 


qn/2 
An = 772 


(n/2)!" 


Our derivation of this formula for a, was a bit clumsy in one sense: it used induction on 
n and produced two formulas depending on the parity of n. The following exercise produces 
a, directly, not relying on induction. It even produces the value for I'(1/2) as a by-product. 
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PROBLEM 10-53. Consider the n-dimensional Gaussian function 


2 2 
—2{ 7-27, 


JJ anf] 2 
ele? = ¢ 


In the following calculations do not worry about the improperness of any of the integrals. 
(This is in accord with the “conservation of mass” principle given on p. 9-25.) 


a. Use Problems 10-49 and 52 to show that 


i el"dx = na, (n/2)/2. 


b. Use Fubini’s theorem to show that 


i ele" de = IL / e tk dap. 
nm k=1 —co 


c. Combine the preceding formulas to achieve 
n 
Gin = (r(a/2))" /v (= } 1) 


d. Use the case n = 2 to show that ['(1/2) = \/z, and thus to obtain the desired 
formula for an. 


THE BETA FUNCTION. There is another important special function, which is closely 
related to the gamma function. These functions deserve attention because they appear so 
often in applications. Again, the word “beta” has no special significance in this context. 
Incidentally, the letter B in the definition is to be regarded as an upper case beta. 











DEFINITION. The beta function is the function (0,00) x (0,00) —> R given by the 
expression 





1 
B(a,b) = | med (ne ai 
0 
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PROBLEM 10-54. Prove that the defining integral for B(a, b) is finite if a > 0 and 
b> 0. Prove that B(b,a) = B(a,b). Also prove that 


a/2 
Bia,0) = 2 | sin?*~! @ cos”?! 6d0. 
0 


The next exercise gives a truly wonderful calculation of the beta function in terms of the 
gamma function. 


PROBLEM 10-55. Explain the validity of each step in the following calculation (you 
may ignore the problems connected with the convergence of the integrals). For any a > 0, 


b> 0, 
| ge ds - | fea 
0 0 


521 po-1e (54) det 


/ 

= | / (eat es deat 
o Jt 
I 


T'(a)P(b) 


(a — £)0 10 16 dadt 


Thus we have found the formula 





['(a)P(b) 


B(a,b) = Tats) 











Notice from Problem 10-54 the special value B(1/2,1/2) = 7. Setting a = b = 1/2 in the 
formula we have just derived gives therefore [(1/2) = \/z. This marks our third derivation of 
the Gaussian integral. 


The formula for B enables us to compute many difficult integrals very quickly. It doesn’t 
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help in computing indefinite integrals such as 


/ sin!” 6d; 


these are sometimes called incomplete beta functions. But the complete beta function is merely 


m/2 
; sin’? 6d6 
0 


5 B(13/2, 1/2) 


T(13/2)P (1/2) 
or (7) 
wget. : 5a rt 


5.3 
2°2 
2-6! 











11-7-37 
pil . 


PROBLEM 10-56. Use the trigonometric version of B(a,a) (Problem 10-54) and the 
formula sin 20 = 2sin @cos@ to derive the duplication formula of Legendre: 


g2a-1 


Pa) = ST (« 4 5) 





PROBLEM 10-57. Use the preceding problem to show that 
1 a 
ri-|JP{-) = 2 


There is another interesting formula for the gamma function, which we can now derive 
rather easily. It is based on the familiar representation of the exponential function, 


se ie (1+=)". 
n 


nC 


Now consider the integral representation of (a), but instead of e~ in the integrand insert the 
polynomial (1 — t/n)”". The resulting integral that we consider is 


n t n 
rT, = | as (1-=) dt. 
0 nm 
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Replace t by ns to get 
1 
= nf s* 1 (1 — 8)" ds. 
0 


This integral can be evaluated explicitly simply by integration by parts n times. We can 
actually circumvent that by noting 


Tr, = n°B(a,n+1) 
pohalin +) 





T(a+tn+1) 
_ yt T(a)n! 
~~ (a+n)---(a+ lal (a) 
nn! 





(a+n)---(a+l)a 
It is not difficult to show that [,, has the limit (a) as n — oo. Thus we obtain 
THEOREM. For0<a<o, 





nn! 
Tr = ii ; 
(2) care (a+n)---(a+l1)a 














To be completely honest, we need to provide a justification of the equation 
hin Ty, = Fla): 


There are several ways to prove this, and here is an ad hoc one that has the advantage of 
being elementary. First note that 


e” > 1-2 for all O< a <1. 


Thus 


Therefore we immediately see that 


T'(a) > / oe “a 
0 


V 
i=) 
3 
oH 
a 
iF 
aN 
_ 
| 
Silo 
NY 
3 
Q 
oe 
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Thus we need to establish some sort of reverse inequality. 
One such inequality is this: 


g(t) =e (a- ~\ "4 =| ; 


We need to show that g(t) > 1. Since g(0) = 1, it will suffice to prove that g’(t) > 0. But 


EA . & 
ge =< (1-=) +—- for O<t<n. 
n n 


To prove this define 








n= e[(-t)+t-(0-8) 
2s Eft (0-2) +2 
> 0. 


Now let € > 0 be arbitrary, and choose to sufficiently large that 


/ foe df =e 
to 


to 
e+ f fled 
0 
to a to t 
c+ f p> (1-=) a+ | (ar 
0 n 0 n 


e+ [ ao Oe 
0 n (a+1)n 


a+1 
to 


(a+1)n 


We conclude that for all n > to 


IA 


['(a) 


IA 


IA 


= e+ dl aye 


Therefore, for all sufficiently large n, 
Tia) < 2e4+T,. 


We conclude that I’, has the limiting value [(a). 
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PROBLEM 10-58. Use the above theorem to prove 








= 2 24 a 2n 2n 
Wallis’ formula 5 Pet 133 5 Qn -1 2n+1 











(This is usually written 





=~ 
cone 
eo | & 
on} aS 
CUS 
aes 


I. Notation for the Jacobian determinant 

There is a very fine notational device which aids the memory in handling changes of 
variables. Returning to the theorem of Section G, where y = ®(x) denoted the variable 
change, we have the scale factor 


det ®'(x). 
The clever notation for this quantity is 
ON, ere ,Yn) 
O(x1, asia’ sie) 


That is, 

Oy, nays Un) 

—— = = det(0y;/0z;). 

.a)° ee 
With this notation the formula for changing variables looks almost like that of the one- 
dimensional case (except for the absolute value sign): 


O(y1, ree Vin) 


O Digaess ty a 


fay = [  (Q)) | 


ye Ag re Ay 





PROBLEM 10-59. Show that 


Oi ens) iy OMiesan4 the) 
O(Y1,---5 Yn) O 


EXAMPLE. Let A be the region in the first quadrant of the x — y plane determined by the 
inequalities x < y < 3” and 2 < xy < 4. For this region we would likely be interested in a 
variable transformation of the sort 
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u =y/a, 
Uv = xy. A 


> 


Then 











O(u, v) ae es ) 


y x 
—2y/ax. 


Notice that we can express x, y explicitly in terms of the new variables as 


G = f=, 
‘; = Juv. 
Therefore the change of variables we obtain is 
dady = ~ dudv 
2y 


1 
= 5 tude 


and we obtain for example 


| sdzdy = / aids 
A u 2U 


l<u<3 
2<vu<4 


1 3 4 
aS / | Jou ?/*dudu 
1 J2 


1 3 4 
= a wPau- f Judu 
21 2 


2 (1- a) 38 2v2) 


(1- =) (4— V2). 


wl Ee Dl FR 


PROBLEM 10-60. Calculate the area of the set A in the preceding problem. 


(Answer: log 3) 
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PROBLEM 10-61. For 2 x 2 matrices A regard det A as a function from R* to R, 
and then calculate 


a. the average of | det A] over the cube (0, 1]*; 


b. the average of | det A] over the cube [—1, 1]*. 


(Answers: 13/54 and 10/27) 


PROBLEM 10-62. Repeat the preceding problem, only with | det A| replaced with 
(det A)?. 


(Answers: 7/72 and 2/9) 


PROBLEM 10-63. For 3 x 3 matrices A regard det A as a function from R® to R, 
and then calculate 


a. the average of (det A)? over the cube (0, 1]°; 


b. the average of (det A)? over the cube [—1, 1]°. 





(Answers: 5/144 and 1/9) 


My student Stefan Allan recently showed me a problem from page 187 of Niven, Zuckerman, 
and Montgomery, An Introduction to the Theory of Numbers, 5th edition, 1991, and the 
following problem is based on that one. 


PROBLEM 10-64. The floor function defined on R is defined by the expression 
|t| = largest integer less than or equal to t. 


(Sometimes this is designated [t].) Let a € R” have integer coordinates, and evaluate the 


| |aex|dx. 
[0,1]" 


(HINT: change variables with x; = 1 — y;.) 


integral 
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PROBLEM 10-65. Under the same conditions as in the preceding problem, show 
that 


/ laex|dx = —1. 
[-1,1]" 





